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I.  THE  DEMONSTRATION 


Ob  j ec  t i ve 


Obtain  higher  resolution  in  the  doppler  frequency  domain  for  improved 
target  detectability  in  distributed  clutter. 

Approach 

Intra-pulse  frequency  diversity  waveform  is  used  to  increase  the  number  of 
independent  data  samples  for  pulse-doppler  processing  without  increasing  the 
observation  length.  These  independent  data  samples  are  coherently  processed 
via  the  newly  developed  C ' '  3  M  F  B  L  P  method  [1]  to  achieve  a  stable  high- 
resolution  doppler  spectrum. 

Data  Acquisition  and  Processing 

The  RADC  L-band  facility  and  AP-120B  array  processor  were  used.  See  (2J 
for  details  of  the  waveform  design,  configuration  and  programs.  Matrix  opera¬ 
tions  were  performed  on  the  complex  baseband  radar  data  to  provide  coherent 
processing  of  the  doppler  spectrum  by  eigenvalue  decomposition. 

Contents  of  This  Demonstration 

(1)  The  intra-pulse  frequency  diversity  waveform  which  drives  the  L-band 

transmitter:  three  subpulses  of  three  microsecond  each  with  the  carrier 

frequencies  of  -6,  -1  and  +4  MHz  (with  respect  to  the  reference  fre¬ 
quency)  . 

(2)  The  spectrum  of  the  above  waveform. 

(3)  The  range-doppler  plot  witli  t lie  conventional  FFT  based  pulse-doppler  pro¬ 
cessing  (8  pulses  with  rectangular  window):  the  weak  target  in  the  Lange- 
extended  strong  clutter  is  not  detectable  since  the  doppler  frequency 
separation  is  smaller  than  the  Fourier  resolution  limit. 

(4)  Same  as  (3)  above  with  Hamming  window  for  lower  sidelobes:  target  is 
still  not  detectable. 

(5)  The  range-doppler  plot  with  the  new  method:  the  target  is  clearly  detec¬ 
table. 

(6)  The  "zoomed-in"  range-doppler  plot  with  the  new  method. 
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A  Comparison  of  Smoothing 
Approaches  for  Angie 
Measurement 

JUN-XIOWC  ZHU,  Student,  IEEE 

HONG  WANG,  Member,  IEEE 
Syracuse  University 


For  Ugh-rt solution  active  direction  finding  of  completely 
correlated  targets,  one  can  use  two  approaches,  Le„  the 
frequency-domain  smoothing  (1)  and  the  non-frequency -domain 
smoothing  such  as  various  spatial  smoothing  (2-7),  or 
multidimensional  search  of  signal-subspace  (8,  9],  A  performance 
congparfcon  of  the  two  approaches  for  the  fluctuating  target 
cases  under  equal  transmitted -energy  constraint  is  presented. 

Both  theoretical  analysts  and  simulations  are  used  to  study  the 
performance  of  the  detection  (determination  of  number  of  targets) 
and  angle  estimation.  It  la  found  that  the  frequency-domain 
smoothing  can  significantly  outperform  the  non- frequency -domain 
smoothing  under  the  equal  Irana  milled -energy  constraint. 


Manuscript  received  November  6.  1987 
IEEE  Log  No,  28658. 

This  work  was  supported  by  Rome  Air  Development 
Center  under  Air  Force  Contracts  F30602-81-C-0169  and 
F30602-88-D-0027(A-8-ll24),  and  by  RADC-LDF 

Authors'  address:  Dep  t,  of  Electrical  and  Computer  Engineering, 
Syracuse  University,  Syracuse,  NY  13244-1240. 


i)0 18-925 1/894)700-05 29  Si.00  ©  1989  IEEE 


High-resolution  direction  finding  of  angularly 
closely  spaced  targets  has  received  increasing  attention 
in  the  past  few  years.  Many  methods  have  been 
suggested  and  are  still  under  further  investigation,  e.g., 
MEM,  MLM,  MUSIC,  modified  FBLP  (Min-Norm), 
ESPRIT,  etc  Unfortunately,  the  detection,  resolution, 
and  estimation  performances  of  these  methods  often 
suffer  severe  degradation  in  the  cases  of  completely 
correlated  target-signal  returns  caused  by  severe 
multipath  or  “smart”  jamming. 

The  degradation  can  stem  from  the  fact  that 
the  covariance  matrix  of  the  target  signals  is  nearly 
singular  or  even  completely  singular.  One  class  of 
methods  to  remove  such  singularity  is  the  so-called 
spatial  smoothing  [2-7j.  For  example,  the  full  sensor 
array  can  be  divided  into  partially  overlapping 
subarrays  to  average  the  covariance  matrices  of  the 
subarray  output  vectors.  It  is  obvious  that  such  a 
spatial  smoothing  method  can  only  be  applied  to 
the  arrays  with  uniformly  spaced  identical  sensors. 
Besides,  a  trade-off  has  to  be  made  for  a  given  array 
between  the  size  of  the  subarrays  and  the  number  of 
subarrays,  since  the  final  performances  depend  on  both 
parameters.  Moreover,  if  the  number  of  sensors  is  not 
large  enough,  the  number  of  s>  .arrays  for  the  spatial 
smoothing  might  be  too  small  to  effectively  remove 
the  singularity.  One  can  also  easily  see  that  such  a 
spatial  smoothing  method  docs  not  make  full  use  of 
the  information  contained  in  the  data  vectors  since  the 
off-diagonal  blocks  of  the  estimated  covariance  matrix 
are  discarded.  Other  methods  to  solve  the  difficulty 
include  the  maximum  likelihood  parameter  estimation 
[10),  which  in  general  requires  a  computationally 
extensive  multidimensional  search  algorithm  except 
in  the  case  of  the  symmetric  multipath,  and  the 
multidimensional  search  of  signal-subspace  [8,  9], 
which  is  also  computationally  demanding. 

For  active  direction  finding  systems  in  which  one 
has  control  of  the  transmitted  waveform,  Wang,  Li, 
and  Zhu  [1]  proposed  the  frequency-domain  smoothing 
approach  in  contrast  to  the  above  mentioned  methods 
which  we  refer  to  as  non-frequcncy-domain  smoothing. 
The  frequency-domain  smoothing  approach  is  based 
on  the  use  of  frequency  diversity  signaling  and 
the  coherent  processing  idea  which  was  originally 
developed  in  (11-13)  for  passive  direction  finding  of 
wideband  sources. 

We  present  a  performance  comparison  of  the 
two  approaches  for  active  direction  finding  of 
completely  correlated  fluctuating  targets  using 
multiple  observations  (snapshots).  A  frequency 
diversity  signaling  with  the  equal  transmitted -energy 
constraint  is  assumed  for  the  frequency-domain 
smoothing  in  order  to  have  a  fair  comparison. 

This  paper  is  organized  as  follows.  In  Section  II 
we  set  up  the  data-vector  models  for  fluctuating 
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Fig.  1.  Frequency  diversity  signaling  pulse. 

targets  and  introduce  our  notations.  In  Section  III 
we  compare  the  potential  performances  of  angle 
estimation  of  the  frequency-domain  smoothing  and 
non-frequency-domain  smoothing,  via  evaluating  the 
Cramdr-Rao  Lower  Bounds  (CRLBs)  associated  with 
the  two  approaches  under  the  equal  transmitted-energy 
constraint.  Such  a  study  is  to  show  the  effectiveness 
of  the  two  smoothing  approaches  with  regard  to  the 
estimation  performance  in  an  algorithm-independent 
way  for  generality  of  the  conclusion.  We  then  compare 
in  Section  IV  the  estimation  performance  of  a  typical 
frequency-domain  smoothing  algorithm  with  that  of  a 
typical  non-frequency-domain  smoothing  algorithm, 
using  statistical  simulations.  Section  V  is  devoted  to 
the  detection  performance  comparison,  followed  by  a 
summary  of  conclusions  and  discussions  in  Section  VI. 


II.  DATA  SET  MODELING 

Two  signaling  schemes  are  used  in  this  paper. 

The  frequency-domain  smoothing  approach 
uses  a  frequency  diversity  signaling,  while  the 
non-frequency-domain  smoothing  approach  uses  a 
conventional  (no  diversity)  signaling. 

A.  Frequency  Diversity  Signaling 

We  consider  the  following  simple  frequency 
diversity  signaling  as  illustrated  in  Fig.  1,  where  each 
subpulse  of  width  r'  =  t/J  has  a  different  carrier 
frequency  with  an  arbitrary  narrowband  modulation  of 
its  own  For  convenience  of  presentation  we  assume 
that  the  carrier  frequencies  fr  j  =  1,2,...,/  are 
uniformly  spaced,  i.e.,  /j+i  -  fj  -  A/  and  that  ail 
subpulses  have  the  same  modulation. 

We  consider  a  linear  array  of  M  wideband  sensors 
covering  the  frequency  band  over  which  the  frequency 
diversity  signaling  is  used.  A  sequence  of  N  identical 
pulses,  each  of  which  contains  /  subpulses  as  in  Fig. 

1,  is  transmitted  and  reflected  by  fluctuating  targets. 
Assuming  that  the  d  targets  in  the  same  range  cell  and 
with  the  same  Doppler  frequency  shift  have  little  range 
change  over  a  time  interval  of  (N  -  l)/pulse -re  petition 
frequency  (PRF),  and  that  /  matched  filters  are 
used  at  all  M  sensors  for  preprocessing,  we  have 
for  each  of  the  N  pulses,  a  set  of  /  M  x  1  data 
vectors  for  processing  over  this  range  cell  in  the 
form  of 

x,  =  A,s,  +  w;,  j  =  1,2 . /  (l) 


where  A.,  M  *  d.  is  the  directum  matrix  at  in¬ 
frequency  fr  sr  d  x  1,  the  target  complex  amplitude 
vector  at  fj  and  w;,  M  x  1.  the  noise  component  from 
the  M  sensor  receivers  of  the  y th  subband  [1], 

We  assume  that  the  frequency  separation  A /  is 
chosen  to  be  large  enough  such  that  sy,  j  =  1,2,...,/ 
are  statistically  independent  [14],  and  that  w;,  j  = 
1,2,...,/  have  zero  mean  and  the  same  correlation 
matrix  <r2I,  where  I  denotes  the  M  x  M  identity  matrix. 
Both  w j  and  s;  are  modeled  as  complex  normal  vectors 
in  this  paper.  The  target  fluctuation  is  assumed  to  lead 
to  the  independence  of  the  N  samples  of  xt  which  are 
identically  distributed.  Thus,  the  data  set  for  processing 
consists  of  N  independent  identically  distributed 
(IID)  samples  of  Xj,  j  =  1,2,...,/,  i.e.,  a  total  of 

iV  x  /  M  x  1  vectors.  We  note  that  x;,  j  =  1,2 . / 

are  not  statistically  identically  distributed  since  A,, 
j  =  1,2,...,/  depend  on  the  carrier  frequencies  fr 
j  =  1,2,...,/. 

Define  the  signal- to- noise  ratio  (SNR)  at  the  / th 
subband  for  the  Icth  target  as 

SNRy(fc)  =  £{|r;(fc)|2}/var{w;(m)}, 

j  -  1,2 . /,  k  =  1,2,. ...d.  (2) 

For  convenience  of  discussion  we  assume  here  that 
the  mean-square  reflection  power  of  the  target  does 
not  depend  on  the  carrier  frequency,  i.e.,  fflryfc)!2} 
independent  of  j.  Since  var{  wy(m)}  =  <r2,  we  have 

SNR,(*)  *  SNR(fc),  j  =  1,2 . /.  (3) 

We  concentrate  on  the  cases  of  two  completely 
correlated  targets  (d  =  2)  with  different  angles  of 
arrival  (AOAs)  8\  j*  0?  and  the  second  arriving  signal 
being  a  delayed  version  of  the  first,  i.e., 

s;(0  =ps\(t-t0)  (4) 

where  p  is  a  real  constant  and  r0  the  path  delay  much 
smaller  than  that  of  the  range  cell.  Let  <?,(1)  and  o,( 2), 
j  =  1,2,...,/  be  the  phases  of  the  target  complex 
amplitudes  s,  ~  [-5,(1)  (2)]T  at  fr  From  (4)  it  is  seen 

that 

o;(2)  =o,(l)-2  7  =  1,2. . /  (5) 

i.e.,  o;(l)  and  o;(2)  are  linked  by  a  constant,  in 
contrast  to  the  two  independent  target  reflection  cases 
where  ©;(1)  and  0:( 2)  can  be  modeled  as  independent 
random  variables. 

For  convenience  of  discussion  we  assume  that 
the  sensors  arc  identical  and  uniformly  spaced  with 
separation  D  —  Ao/2  where  A0  is  the  wavelength 
corresponding  to  /o,  the  central  frequency  of  the 
whole  frequency  band.  The  array  reference  point  is 
chosen  to  be  the  center  of  the  aperture  with  all  AOAs 
referenced  to  the  broadside  of  the  linear  array.  The 
Rayleigh  angular  resolution  limit  of  such  an  array  is 
about  Q  =  2/M  rad. 
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Fig.  2.  CRLBs  (101og],)(l/ir5))  with  frequency  diversity 
completely  correlated  targets. 

B.  Conventional  Signaling 

A  conventional  (no  diversity)  signaling  is  assumed 
for  the  non-frequency-domain  smoothing  approach, 
which  has  a  single  carrier  frequency  /o  for  the  whole 
pulse  with  an  arbitrary  narrowband  modulation, 
t.e.,  7  =  1.  Therefore,  the  data  set  for  processing 
consists  of  rV  M  x  I  [ID  vectors.  Imposing  the  equal 
transmitted -energy  constraint  on  both  signaling 
schemes,  we  have  the  SNR  with  conventional  signaling 


For  a  real,  jointly  normal  distributed  data  set,  the 
Fisher  information  matrix  for  bound  calculation  can  be 
expressed  as  [  16) 


[J(«)]u.*  =  jtr 


9m(0) 


d»u  j 


dm  (9) 


d9v 


where  R(0)  is  the  covariance  matrix  as  a  function  of 
the  parameter  set  vector  9,  m(0)  the  mean  vector  as 
a  function  of  9,  and  9U,  9y  denote  the  uth  and  vth 
parameter,  respectively.  For  a  jointly  complex-normal 
distributed  data  set,  the  above  equation  can  be 
modified  as 


[JWkv  =  tr 

+  2Re 


dm(0) 


d6u 


R -'(*) 


dm(9) 

~d6~ 


(7) 


with  Re()  being  the  real  part  of  a  complex  number. 
In  the  case  of  two  completely  correlated  targets, 
the  unknown  parameters  of  the  data  set  with  a 
conventional  signaling  are  the  two  spatial  frequencies 


ujk  =  2ir(D/\0)sin9k,  k  =  1,2  (8) 


SNRo(fc)  =  IOIog/  +  SNR(*),  (dB), 

k  =  1,2,. ...d  (6) 

i.e.,  101og7  higher  than  that  with  frequency  diversity 
signaling.  We  assume  that  all  other  conditions  and 
assumptions  on  targets  and  array  remain  the  same 
for  the  performance  comparison  of  the  two  smoothing 
approaches. 

III.  ESTIMATION  PERFORMANCE  BOUND 
COMPARISON 

Since  there  are  many  ways  to  implement  smoothing 
with  either  approach,  a  comparison  in  terms  of  some 
performance  bound  is  necessary  for  the  generality 
of  the  conclusions  on  the  potentials  of  the  two 
approaches.  A  well-known  estimation  performance 
bound  is  the  CRLB  which  is  relatively  easy  to  evaluate. 
It  is  well  known  that  CRLB  bounds  the  unbiased 
estimators  only.  Though  a  modified  CRLB  for  biased 
estimators  is  available  [15,  ch.  4],  it  is  in  general 
a  very  difficult  task  to  evaluate  the  required  bias 
function  of  the  specific  estimator.  Since  no  available 
method,  to  the  knowledge  of  the  authors,  has  been 
reported  to  have  a  mean-square-error  (MSE)  below 
the  corresponding  CRLB  and  since  the  CRLB  can  be 
fairly  tight  for  SNRs  above  a  certain  threshold  even 
if  the  data  set  is  small,  the  use  of  the  CRLB  for  the 
purpose  mentioned  above  is  justified. 


the  signal  power  ff] ,  and  <rj2,  the  noise  power  a1,  and 
the  phase  difference  between  the  two  signals  which 
depends  on  the  delay  Iq.  The  unknown  parameters 
with  frequency  diversity  signaling  include  the  signal 
powers  o-;,(y),crj2(;),  noise  power  <rz(j),  j  =  1.2,.. .,7. 
the  delay  t0,  and  two  target  AOA  related  parameters. 
We  still  use  u-t  and  -j;  for  these  two  parameters 
since  the  spatial  frequencies  at  fr  j  =  1,2,. ..,7  are 
deterministic  functions  of  these  two  even  if  /0  is  not 
necessarily  equal  to  one  of  /),  j  =.  1,2 . 7. 

Let  az  be  the  CRLB  to  the  variance  of  an  »t 
or  w’2  estimate.  In  the  remainder  of  this  paper  and 
Figs.  2-5,  the  abbreviation  CRLB  is  used  to  denote 
10l°g1o(l/<T5)  as  a  rather  common  practice. 

We  consider  the  situations  where  p  =  1  and  t0  is 
corresponding  to  a  small  fraction  of  the  range  cell,  i.e.. 
fo<  where  denotes  the  largest  bandwidth 

of  the  subpulses.  In  practice,  these  situations  are 
found  to  be  the  most  important  and  difficult  to 
handle  [10,  ch.  3  and  4).  Fig.  2  shows  the  CRLBs 
as  a  function  of  tofo  for  7  =  7  subbands  with  the 
relative  bandwidth  B/f0  =  10  percent  and  30  percent, 
respectively.  The  CRLBs  oscillate  with  a  damping 
factor  proportional  to  B//0  but  settle  at  a  level 
independent  of  B/f0.  In  contrast,  the  dotted  lines 
in  Fig.  2  indicate  the  maximum  and  minimum  of  the 
CRLB  for  the  non-frequency-domain  smoothing,  which 
would  oscillate  in  between  with  the  period  equal  to 
0.5.  From  this  figure,  we  can  see  that  the  CRLB  for 
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Fig.  3.  CRLBs  (10loglo(l/<r£))  with  frequency  diversity  versus 
number  of  subbands  7:  completely  correlated  targets,  Bf f o  *  10 
percent. 
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Fig  4  CTES*  (!<?!cg10'l/.r' ''  •*"!h  frequency  divr.jic,-  versus 
number  o(  subbands  /:  completely  correlated  targets.  Bj  fa  =  30 
percent. 

the  frequency-domain  smoothing  can  be  independent 
of  rnf, y  as  long  as  enough  bandwidth  is  used.  For  the 
example  of  B/f0  =  30  percent,  when  t0ft>  >  2.5  the 
hound  stays  at  a  level  close  to  the  best  value  of  the 
bound  for  the  non-frcquency-domain  smoothing.  And 
for  fi//o  =  10  percent,  the  bound  gets  into  the  stable 
region  for  r0/o  >  7.5. 

The  effects  of  the  number  of  subbands  J  on  the 
CRLB  for  the  frequency-domain  smoothing  approach 
can  be  seen  from  Figs.  3  and  4  for  the  two  completely 
correlated  targets  with  B/ fo  =  10  percent  and  30 
percent,  respectively.  With  small  J  the  CRLBs  still 
depend  on  the  values  of  r0/o •  When  /  is  large  enough, 
however,  the  CRLBs  can  settle  at  a  level  independent 
of  tofo.  Moreover,  a  slightly  larger  /  is  needed  with 
larger  relative  bandwidth  B/ fo  for  the  CRLB  to 
become  independent  of  tofo- 

From  the  above  study  of  CRLBs  for  the  two 
approaches  we  can  conclude  the  following. 


1)  The  potential  AOA  estimation  performance  ot 
the  frequency-domain  smoothing  approach  is 
expected  to  always  be  better  than  that  of  the 
non-frequency-domain  smoothing  approach. 

2)  The  larger  the  bandwidth,  the  better  the  potential 
AOA  estimation  performance  can  be  of  the 
frequency-domain  smoothing  approach. 

3)  Only  a  small  number  of  subbands  J  are  necessary 
for  the  potential  AOA  estimation  performance  of 
the  frequency-domain  smoothing  approach  to  be 
independent  of  the  unknown  delay  r0. 

In  the  next  section  we  study  the  estimation 
performances  of  a  typical  frequency-domain  smoothing 
method  and  a  typical  non-frequency-domain  smoothing 
method. 

IV.  ESTIMATION  PERFORMANCE  OF  TWO 
SMOOTHINC  METHODS 

For  simplicity  we  choose  the  MUSIC-based 
coherent  signal-subspace  processing  method  described 
in  (11,  17]  and  the  MUSIC-based  spatial  smoothing 
method  of  [4],  Since  there  is  little  analytical  result 
available  atout  the  estimation  performance  of  the 
spatial  smoothing  method,  we  conduct  the  comparison 
via  statistical  simulation.  The  following  parameters  are 
set  up  for  the  simulation. 

1)  the  number  of  targets  d  =  2  with  p  =  1,  i.e., 

SNR(l)  »  SNR(2); 

2)  the  relative  bandwidth  B/f0  =  30  percent  with  the 
number  of  subbands  J  =  7; 

3)  the  number  of  sensors  M  =  8; 

4)  the  spatial  frequencies  u’i  =  0.8227  and  wj  =  0.430 
which  correspond  to  an  angle  separation  equal  to 
0.5fi; 

5)  the  subarray  size  L  =  6  for  the  spatial  smoothing 
method  in  (4); 

6)  the  required  preliminary  estimate  of  approximate 
center  of  the  spatial  frequencies  vith  'he 
frequency-domain  smoothing  method  of  [11]  J  = 

(-T  +w’2)/2; 

7)  the  equal  transmitted-energy  constraint  as  of  (6); 

8)  the  number  of  pulses  iV  =  32. 

Fig.  5  compares  the  mean  square  errors  (MSEs) 
of  the  two  smoothing  methods  for  at  various 
SNR  with  tofo  =  6,  6.25,  and  6.5,  respectively. 

Again,  101og10(l/MSE)  is  actually  plotted.  Fifty 
independent  runs  are  used  at  each  SNR.  The  MSEs 
of  frequency-domain  smoothing  with  the  three  different 
values  of  tofo  are  found  to  be  approximately  the  same 
and  thus  only  the  one  with  tofo  =6  is  plotted  in  Fig.  5. 
Also  included  in  Fig.  5  are  the  corresponding  CRLBs. 
The  frequency-domain  smoothing  method  is  seen  to 
significantly  outperform  the  spatial  smoothing  method. 
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Fig  5  MSEi  of  vi,  with  frequency -domain  smoothing  (B/ ft,  = 
30  percent,  }  m  7)  and  spatial  smoothing  versus  SNR:  completely 
correlated  targets,  M  =*  8,  d  »  2.  -v,  «  0.8227,  *  0  430.  L  »  6. 

iV  -  32. 

and  reaches  the  CRLBs  quite  closely.  We  note  that 
with  ro/o  *  6.25  the  spatial  smoothing  method  fails 
to  reach  the  corresponding  CRLB  even  at  very  high 
SNRs.  A  possible  cause  of  this  may  stem  from  the 
fact  that  with  this  spatial  smoothing  method,  not  all  of 
the  off-diagonal  elements  of  the  estimated  correlation 
matrix  can  be  made  use  of. 

Since  the  frequency-domain  smoothing  method 
needs  a  preliminary  estimate  j3,  it  is  interesting  to 
know  its  effect  on  the  MSE.  It  has  been  shown  in 
[17,  18]  that  the  performance  of  the  frequency-domain 
smoothing  approach  is  not  sensitive  to  this  preliminary 
estimate. 

V.  DETECTION  PERFORMANCE  COMPARISON 

The  above  two  MUSIC  based  smoothing  methods 
both  need  to  first  determine  the  number  of  targets 
d  before  performing  AOA  or  spatial  frequency 
estimation.  One  may  use  AIC  or  MDL  types  of 
methods  [19,  20]  to  conveniently  accomplish  this. 

For  completely  correlated  targets,  (21)  gives  a  spatial 
smoothing  based  modification  of  the  A1C/MDL 
method,  while  a  frequency-domain  smoothing  based 
modification  is  presented  and  analyzed  in  [17],  In 
this  section  we  compare  the  probability  of  correct 
determination  of  the  number  of  targets  with  the  two 
modified  methods  for  two  completely  correlated, 
closely  spaced  fluctuating  targets  under  the  constraint 
of  equal  transmitted-energy.  An  analytical  detection 
performance  study  can  be  found  in  [17]  for  the 
frequency-domain  smoothing  based  method,  but  a 
corresponding  study  for  the  spatial  smoothing  based 
method  is  not  available  to  our  knowledge.  Therefore, 
we  conduct  the  detection  performance  comparison  via 
simulation.  The  parameters  used  in  this  simulation 
are  the  same  as  those  in  the  last  section.  Since  the 
probability  of  underestimating  the  number  of  targets 
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Fig.  6.  Detection  performances  of  spaual  smoothing  and 

frequency-domain  smoothing  methods,  completely  correlated 
targets,  M  -  8.  <1  =  2,  -i  «  0  822'.  =  0  430.  L  »  t>.  S  =  32 

is  the  dominant  detection  error  for  the  cases  of 
completely  correlated  and/or  closely  spaced  signals, 
and  since  our  previous  experience  with  detection 
performance  analysis  has  resulted  in  favor  of  choosing 
the  AIC  rather  than  MDL  penalty  function  (17,  22,  23), 
we  use  the  AIC  penalty  function  here. 

Fig.  6  shows  the  probability  of  correctly 
determining  the  number  of  targets  with  the  spatial 
smoothing  based  method  (21)  and  frequency-domain 
smoothing  based  method  [17]  under  equal 
transmitted-energy  constraint,  using  50  independent 
runs.  In  this  figure,  the  probabilities  for  the  spatial 
smoothing  based  method  with  t0fo  =  6,  6.25.  and 
6.5,  respectively,  are  plotted  as  a  function  of  SNR. 

For  the  frequency-domain  smoothing  based  method, 
the  corresponding  three  curves  are  almost  the 
same.  So  only  the  one  with  tn  f0  =  6  is  included.  The 
detection  performance  of  the  spatial  smoothing  based 
method  is  very  sensitive  to  the  values  of  fo/o  (or  the 
phase  difference  0  =  2xf0/o).  Though  the  detection 
peiformancc  of  the  frequency  d'''r?in  based  method 
is  measured  about  2  dB  poorer  in  terms  of  SNR  than 
the  spatial  smoothing  based  method  with  i0/o  =  6.5,  it 
performs  significantly  better  than  the  other  two  cases 
of  the  spatiai  smoothing  based  method.  Therefore,  one 
may  conclude  that  the  frequency-domain  smoothing 
based  detection  method  is  more  able  to  provide  a 
well-performed  solution  for  determination  of  the 
number  of  completely  correlated,  closely  spaced 
fluctuating  targets. 

VI.  CONCLUSION  AND  DISCUSSION 

For  high-resolution  active  direction  finding  of 
completely  correlated  fluctuating  targets  the  results 
are  in  favor  of  using  the  frequency-domain  smoothing 
approach,  which  consists  of  the  frequency  diversity 
signaling  and  a  coherent  wideband  processing 
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algorithm  such  as  the  one  described  in  ill.  l’|, 
rather  than  the  nor.-frequency-domatn  smoothing 
approach  wht>  ■  uses  the  conventional  (no  diversity) 
signaling  ■  .  a  narrowband  processing  algorithm 
such  is  the  spatial  smoothing  method  of  [4,  21].  The 
study  of  the  estimation  performance  bounds  also 
indicate  that  the  frequency-domain  smoothing  can 
be  expected  to  achieve  a  much  better  performance 
than  any  non-frequency-domain  smoothing,  especially 
for  an  array  with  a  small  number  of  sensors.  We 
note  that  for  nonfluctuating  targets,  (24}  contains  a 
study  which  leads  to  the  same  conclusion  in  favor  of 
using  the  frequency-domain  smoothing.  We  also  note 
that  the  frequency-domain  smoothing  can  be  used 
with  atTays  of  non-uniformly-spaced,  non-identical 
sensors.  It  should  be  pointed  out  that  passive  direction 
finding  of  completely  correlated  wideband  sources 
has  been  successfully  solved  by  using  the  idea  of  the 
frequency-domain  smoothing  (11,  17]. 
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Smite  remarks  on  the  above  analysis  are  in  order  now 

Ht  nulls  !  .Since  the  estimate  Plui,.  u><).  gi veil  by  1 4. > .  (K).  and 
i'll,  satislies  eondition  I )  lor  t  /.  k  I  6  S!  and  condition  2)  lor  I  k  |  > 
A.  respectively,  we  can  conclude  that  the  estimate  Pi  w,.  w_, )  given 
by  Kiniura  and  Honoki  s  hybrid  approach  coincides  with  the  true 
ME  estimate  when  [l„t  ca, )  is  positive  dehniic.  and  the  elements  of 
the  lirst  row  ol  its  inverse  fl„  1  ( iu,  )  satisfy  (26). 

Remark  2  If  the  original  estimate  II  (  u>,  )  is  of  Tocplitz.  H  (<j| ) 
is  directly  taken  as  II, I  w,  );  in  the  other  case.  II„(  w,  )  is  obtained 
bv  averaging  FI  t  u) ,  ) .  see  |l.  oq.  <  20  j  | .  In  a  way  similar  lo  our 
above  analysis,  it  can  he  straightforwardly  concluded  that  for  a 
multichannel  process  .Vt  t ).  the  power  spectrum  estimate  [I„l  ca,  I  is 
the  (rue  ME  estimate  when  every  element  |  II  ,  1  ( w,  )  |1(  of  the  in¬ 
verse  n,/‘(u),  )  has  the  lorm  ol  12, F  C',‘  11  exp  (  yea,/  )  lor  1/  |  > 
L  This  can  be  regarded  as  another  font)  of  (26)  for  (he  multichan¬ 
nel  1-D  process  II  is  worthwhile  to  point  out  even  if  the  estimate 
II.,  I  w,  1  is  the  true  ME  estimate,  (he  linal  estimate  Pi  u>,.  ca.  I  ob¬ 
tained  by  using  II  , ( ca,  )  is  not  necessarily  the  ME  estimate,  as  will 
be  stated  in  the  next  remark 

Remark  J:  Since  the  positivity  of  Fl,,  ( ca , )  is  available  for  the 
ease  ol  cyclic  and  skew-cyclic  Tocplil/  Ri  I  ).  Kimura  anil  Honoki 
have  cniiiccliircd  that  their  linal  estimate  Pi  w,.  w. )  coincides  with 
the  true  ME  estimate  lor  such  a  case  The  above  analysis  conclu¬ 
sion  shows  that  Kimura  and  Honoki's  conjecture  is  not  true  in  gen¬ 
eral  since  the  cyclic  and  skew-cyclic  Tocplitz  Ri  I  )  cannot,  in  gen¬ 
eral.  guarantee  that  all  the  elements  of  the  lirst  row  ol  the  inverse 
l!„  'tui,  )  satisly  (26).  To  show  this  point,  we  present  a  counter¬ 
example  ol  Kimura  and  Honoki's  conjecture. 

It  is  always  possible  to  lind  a  multichannel  process  ,V(r)  such 
that  it  has  the  following  power  spectrum: 


is  tlie  sum  ol  mlinite  terms;  this  implies  that 

i"  i"  exp  j  jiui,l  +  u)>k  )  | 

I  1  a -  i/ui,  ,/u--  at  t). 

J  rJ-r  r  (  W, ,  uli  ) 

for  k  <  K.  |  / 1  >  L. 

IV.  Conclusions 

Some  important  theoretical  problems  remain  to  be  solved  in  Ki¬ 
mura  and  Honoki's  hybrid  approach  to  high-rcsolutu.n  2-D  spec¬ 
trum  analysis.  The  most  important  problem  is  probably  to  know 
when  the  linal  estimate  is  close  to  the  true  ME  estimate.  In  this 
correspondence,  we  have  analyzed  such  a  problem,  and  our  wor) 
has  provided  a  "theoretical"  and  “practical  "  solution  to  this  prob¬ 
lem.  Using  our  result,  one  can  easily  know  if  the  linal  estimate 
given  by  Kimura  and  Honoki's  hybrid  approach  does  coincide  with 
the  true  ME  estimate. 

We  have  also  presented  a  counterexample,  and  have  shown  that 
Kimura  and  Honoki's  conjecture  on  the  above  problem  is  not  true 
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FKui,) 


_ 5  -t-  2  cos  ui, _  -2  sin  ca, 

25  +  20  cos  ca,  -t-  4  cos  ( 2w, )  25  +  20  cos  w,  +  4  cos  ( 2ca, ) 

_ -2  sin  co | _  _ 5+2  cos  co i _ 

_25  +  20  cos  ca i  +4  cos  (  2ca, )  25  +  20  cos  w,  +  4  cos  (  2ui, )  _ 


/„<uh)  /i;(uti) 
/;,  ( ca, )  /..  ( ca, ) 


Hy  using  deli  ml  ion  (  I  I  a),  it  is  easy  to  show  lit  w,  )  is  of  cyclic 
Tocplil/.  Since  there  are  / , ,  I  w,  )  =  /  t  ca,  )  and  / ,  -  (  ca,  I  =  /  .,  (  ui,  ) 
helwcen  the  elements  ol  III  uj,  ).  wc  see  that 

I  '  1 

/?(/)  =  —  \  n ( u>, )  exp  ( /u>,/ ) 

2*  .)  » 


Adaptive  Beamforming  for  Correlated  Signal  and 
Interference:  A  Frequency  Domain  Smoothing 
Approach 


is  ci  I  so  ol  cyclic  Tocplil*.  Hence.  fl(uj,)  is  directly  taken  as 
II ,(  w,  l,  and 


5  +  2  cos  u>j  2  sin  uj, 

ii  1 1 wd, )  —  n  1 1  u/j )  - 

I  2  sin  uj,  5+2  cos  uj( 


(281 


Note  that  the  given  flfuj,)  is  positive  definite  and  is  simply 
shown  to  he  the  ME  estimate  let'  Remark  2)  and  the  corresponding 
Rif  \  is  of  cyclic  Toeplit*.  hut  the  final  spectrum  estimate 
w-  i  obtained  hy  Kimura  and  Honoki's  hybrid  approach  'rom  such 
Iliuj,  i  floes  not  coincide  with  ihc  true  Ml:  estimate  since  lor  ihc 
•  jsc  ol  k  -  l)  and  i  f  i  >  /.. 


uj,  ) h*,, ( UJ,  1 

"  h,HI(  uJ|  ) 

h,9,i  UJ,  )/>,»,{  w,  )  +  />,n  (  uj,  )/>■*  (  uJj  ) 

ft,*,  A  w,  ) 

4  sm'  uj, 

=  5  +  2  cos  ojj  + - 

5  +  2  cos  uj. 
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Abstract— This  correspondence  succests  a  new  approach  to  wide* 
hand  adaptive  hcamforminu  for  correlated  signal  and  interference. 
(Jnltke  other  approaches,  such  as  spatial  smoothing  the  ness  approach 
solves  the  siunal  cancellation  problem  by  employing*  the  idea  of  fre- 
(juency  domain  smoothing.  Advantaces  of  frequency  domain  smooth- 
inu  over  spatial  smoothing  are  identified  in  this  correspondence.  Pre¬ 
liminary  performance  studies  of  a  simple  frequency  domain  sinoothtnu 
method  show  that  proper  spatial  tillering  can  he  achieved  hy  frequency 
domain  smoothing,  svhether  the  desired  siunal  and  the  inlerlercnce  arc 
correlated  or  not. 
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1.  Introduction 

Signal  cancellation  is  a  problem  of  adaptive  beamlorming  (AUK) 
when  the  desired  signal  I  root  the  look  direction  and  the  interference 
arc  highly  correlated.  One  cure  is  the  so-called  spatial  smoothing 
method  1 1 1.  by  which  the  lull  sensor  array  is  divided  into  partially 
overlapping  suharrays  to  enable  an  average  of  the  covariance  ma¬ 
trices  of  the  subarray  output  vectors.  Such  an  average  was  shown 
to  be  able  to  destroy  the  correlation  between  the  signal  and  inter¬ 
ference.  Unfortunately,  spatial  smoothing  can  only  be  applied  to 
arrays  of  uniformly  spaced,  identical  sensors.  Besides,  a  tradeoff 
has  to  be  made  lor  a  given  array  between  (he  si/e  of  the  suharrays 
and  the  number  of  suharrays.  since  (he  interlerenee  rejection  per¬ 
formance  with  spatial  smoothing  depends  on  both  parameters. 
Moreover,  if  the  number  of  sensors  is  not  large  enough,  the  number 
of  subarrays  available  for  spatial  smoothing  might  be  too  small  to 
lead  to  any  sigmlieant  improvement. 

A  rather  different  approach  for  wide-hand  ABF  is  given  in  |2|. 
which  docs  not  sutler  from  the  signal  cancellation  since  it  employs 
the  "complete  nulling"  criterion  at  every  subband.  Such  a  crite¬ 
rion.  however,  may  .sacrifice  too  much  signal  power,  resulting  in 
an  array  output  signal-to-intcrferenec-plus-noi.sc  ratio  (SINR)  lower 
than  that  potentially  achievable,  especially  when  the  interference 
sources  arc  not  very  strong  and  their  angles  of  arrival  arc  close  to 
that  of  the  desired  signal  |3|. 

In  this  correspondence,  we  present  a  new  approach  to  wide-hand 
ABF  lor  correlated  signal  and  interlerenee.  The  key  idea  is  to  per¬ 
form  a  smoothing  operation  in  the  frequency  domain,  which  was 
originally  presented  in  the  contest  of  high-resolution  direction  find¬ 
ing  of  multiple  widc-band  sources  |4|. 

Fig.  I  shows  a  general  configuration  of  a  frequency  domain  im¬ 
plementation  of  wide-hand  bandpass  ABF.  For  simplicity  of  pre¬ 
sentation.  we  consider  a  linear  array.  The  sensors  do  not  need  to 
be  uniformly  spaced,  nor  do  their  patterns  need  to  he  identical. 
Following  each  sensor  is  a  hank  ol  J  narrow-band  tillers  covering 
the  whole  Irequeney  band  of  interest.  The  bandwidth  ol  each  sub¬ 
band  is  assumed  to  be  a  lew  percent  of  its  central  frequency  /,y  = 
1.2.  ■  ■  •  .  J  so  that  for  each  subband,  the  narrow-band  array  out¬ 
put  representation  is  valid.  The  ABF  weight  vectors  tv,,  y  =1.2. 

■  -  •  .  J  arc  updated  according  to  the  chosen  optimization  criterion 
and  interlerenee  environment.  Several  criteria  all  result  in  (he  form 
of 

tv;  =  Or/*,;'<!,(0,)  y  =  i.  2,  •  •  •  . y  { i ) 

where  «,  is  u  criterion  dependent  constant,  R,,  is  the  correlation 
matrix  of  the  complex  output  vector  x, ( l )  of  the yth  suhhand  array, 
and  a;(  9,  )  is  the yth  suhhand  array  direction  vector  at  the  look  an¬ 
gle  9,.  When  the  interference  is  completely  or  partially  correlated 
to  the  desired  signal  from  0, .  the  signal  cancellation  occurs  and  the 
claimed  optimum  performance  is  completely  lost  due  to  the  viola¬ 
tion  of  the  basic  assumption. 

It  should  be  pointed  out  that  if  the  existence  of  the  correlation 
between  the  signal  and  interference  is  known  in  advance,  then  one 
should  lake  advantage  of  the  signal-correlated  interlerenee  to 
achieve  an  even  higher  SINR.  instead  of  simply  trying  to  suppress 
all  interference  |5|.  However,  such  a  priori  knowledge  is  not  al¬ 
ways  available,  nor  always  easy  to  obtain  from  the  received  data. 
Therefore,  the  need  still  remains  for  finding  h>,.  y  =  I.  2.  *  •  •  .  J 
such  that  the  interference  suppression  performance  is  much  better 
ill. mi  that  using  ( I  i  when  llic  signal  and  interlerenee  are  correlated, 
and  close  to  that  using  (  I )  when  the  signal  and  interlerenee  arc  not 
correlated. 

In  the  following  discussion  we  assume: 

1 )  the  angle  of  arrival  of  the  desired  signal  0,  is  available  and 
used  as  the  array  look  direction  9, ; 

2)  the  interference  sources,  correlated  or  uncorrelaled  with  the 
desired  signal,  occupy  the  same  Irequeney  band  as  the  desired  sig¬ 
nal.  and 

3)  the  receiver  noise  is  uncorrelated  sensor  to  sensor 


Fig,  I .  Configuration  of  Irequeney  decomposed  wide  band  heamlormer 


II.  Frit/ui  ncv  Domain  Smiioiuino 

it  is  easy  to  see  that  the  correlation  matrix  of  the  yth  suhhand 
array  output  has  the  lonn 

-  a,*,*!'  +  °;i  (2) 

where  Ar  M  x  J.  J  <  M  is  the  direction  matrix  associated  with 
the  desired  signal  and  il  -  I  interference  sources;  Rr  d  x  </  is  the 
correlation  matrix  of  the  desired  signal  and  inlcrfcrcnec;  and  o;  / 
is  the  correlation  matrix  of  the  receiver  noise  of  the  yth  suhhand. 
When  there  is  a  complete  correlation  between  the  desired  signal 
and  an  interference  source.  Rr  j  =  1 . 2.  ■  ■  •  .  J  becomes  singular 
and  the  so-called  signal  cancellation  occurs. 

As  pointed  out  in  (4|,  the  frequency -domain-smoothed  correla¬ 
tion  matrix  R  =  Rt  is  nonsinguiar  in  general.  This  lact  leads 
to  the  opportunity  of  employing  a  frequency -domain-smoothed 
“correlation  matrix"  to  replace  k,t  ol  (I).  so  that  the  heamlormer 
can  deliver  a  reasonably  good  spatial  filtering  performance  insen¬ 
sitive  to  whether  the  desired  signal  and  interference  are  (c  mi- 
plctcly)  correlated  or  not. 

Parallel  to  the  direction  finding  problem  using  the  idea  of  fre¬ 
quency  domain  smoothing  |4|.  a  class  of  adaptive  wide-hand  beam- 
formers  with  frequency  domain  smoothing  can  be  developed  In 
dilferent  applications  In  this  correspondence,  however,  we  are  only 
interested  in  presenting  one  example  which  is  suitable  lor  a  number 
of  practical  applications  in  radar,  sonar,  and  spread-spectrum  com¬ 
munication. 

To  implement  the  smoothing  operation  on  Rr  j  =  1.2.  •  •  • 

J.  a  frequency  domain  transformation  represented  by  T,  must  he 
performed  on  the  estimate  of  the  correlation  matrix  R ,,  such  that 

T,A,  =  Atl,  y  =  I.  2.  •  •  •  .  J  (3) 

where  A«  is  the  direction  matrix  of  the  array  at  the  central  fre¬ 
quency  fn  (  /i,  may  be  equal  to  one  of/,  y  =  1 .  2.  •  ■  ■  .  J  ).  Ob¬ 
viously,  there  arc  many  possible  choices  of  F,  which  could  achieve 
(3),  but  all  of  them  would  ncetl  the  angles  of  arrival  of  the  inter¬ 
ference  which  arc  not  available  to  the  heamforming  subsystem. 
Therefore,  approximations  of  T,  ol  (3)  have  to  he  used  in  practice 

One  possible  way  is  to  do  a  preselected  angle  approximation  to 
the  perfect  transformation  matrix  Tr  Let  6 ,  be  the  look  direction 
of  the  array,  and  9,„  *  6, .  m  =  1.2.  ■  •  •  .  M  -  I  be  M  -  I 
different  angles  uniformly  covering  the  whole  angular  domain  of 
concern.  Conceptually  summarized  below  arc  the  steps  used  to  im¬ 
plement  the  Ireqticncy-domain  smoothed  adaptive  heamlormer  with 
a  preselected  approximation  ol  the  transformation. 

1)  Form  R,r  the  estimate  of  the  correlation  matrix  R,,.  from  the 
K  sampled  complex  output  vectors  x,(  i ).  j  -  1 .  2.  ■  •  •  .  J: 

2)  perform  the  frequency  domain  transformation  to  obtain 

/ 

R,  =  X  14) 

where 

T,  =  A„A,  '  I  5  I 
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w  ith 

■V  =  (<!.(«,  I  I  o  («„  ,)|;  Ifi) 

3)  form  the  Ircqucncy-domain-smoothed  "correlation  matrix'' 
estimate  lor  each  subbund 

*„==  f T,  ')".  (7) 

4)  obtain  the  beamforming  weight  vector  for  (he  yth  subbund 

u>.  =  *,;'a,(0, ).  j  =  l.  2.  •  •  •  ,  J  (8) 

We  note  that  significant  reduction  in  computing  for  Tt  and  Tt  ' 
can  he  expected  by  exploiting  the  structure  of  A,.  In  the  following 
section,  we  will  present  our  preliminary  performance  study  of  the 
above  method 

III  Preliminary  Performance  Study 

In  this  section  we  will  show  the  beam  patterns  of  both  frequency 
domain  smoothing  and  spatial  smoothing  For  simplicity,  we  con¬ 
sider  an  array  of  M  =  10  uniformly  spaced,  omnidirectional  sen¬ 
sors  with  the  space  between  sensors  equal  to  one-hall  of  the  wave¬ 
length  corresponding  to  J„.  and  we  assume  that  the  desired  signal, 
interference,  and  receiver  noise  all  have  Hat  spectra  in  the  same 
frequency  band  with  a  relative  bandwidth  B/f„  =  0.30.  The  whole 
band  is  covered  by  J  -  10  subbands.  The  angle  of  arrival  of  the 
desired  signal  is  taken  to  be  d,  =  35°.  Two  interference  sources 
are  present.  The  first  one  is  completely  correlated  with  the  desired 
signal,  being  a  delayed  version  with  delay  t„  =  15/%.  The  second 
interference  source  is  a  strong  jammer  which  is  uncorreluted  with 
the  desired  signal  The  angles  of  arrival  of  the  two  interference 
sources  are  30°  and  50°.  respectively. 

In  Fig.  3,  the  solid  line  is  the  beam  pattern  of  the  frequency 
domain  smoothing  method.  The  two  dashed  lines  arc  the  beam  pat¬ 
terns  of  the  spatial  smoothing  method  with  the  subarray  size  equal 
lo  4  and  b.  respectively,  and  ihp  beam  pattern  without  any  smooth¬ 
ing  is  plotted  as  the  (Jotted  line.  The  protection  of  frequency  do¬ 
main  smoothing  against  both  interference  sources  is  seen  to  be  much 
better  than  that  of  spatial  smoothing. 

Since  the  beam  pattern  of  frequency  domain  smoothing  is  the 
sum  ol  ihc  J  subbaml  beam  patterns,  the  sharp  null  at  the  angle  of 
the  strong  interference  indicates  that  the  J  subband  beam  patterns 
musi  have  nulls  almost  at  the  same  angle  even  though  the  transfor¬ 
mation  used  is  just  an  approximation. 

Comparing  the  two  beam  patterns  of  spatial  smoothing,  one  can 
see  that  the  protection  against  the  correlated  interference  can  be 
gradually  improved  by  reducing  the  size  of  the  subarrays.  Such  a 
reduction  ol  the  subarray  sine  for  mure  spatial  smoothing  not  only 
is  limited  by  the  number  ol  interference  sources,  but  also  degrades 
ihc  nulling  performance  at  the  angle  of  the  strong  interference  due 
to  the  loss  ol  angular  resolution  with  the  reduced  subarray  size  |6|. 
In  contrast,  frequency  domain  smoothing  is  free  from  this  problem 
as  shown  in  Fig  2. 

If  the  interference  exists  only  over  purl  of  the  frequency  band, 
( 7 )  and  (8)  indicate  that  this  frequency  domain  smoothing  method 
will  produce  nulls  in  all  J  subband  beam  patterns  toward  (he  par¬ 
tial-hand  interference.  To  avoid  such  unnecessary  nulls  in  the  in- 
tcrlcrcncc-tree  subbands,  one  may  simply  apply  Steps  1 ) -4 )  to  those 
subbands  containing  interference,  and  use  the  conventional  beam¬ 
forming  method  lor  the  interference- free  subbands  If  the  interfer¬ 
ence  sources  have  a  colored  spectrum,  a  more  sophisticated  method 
ot  frequency  domain  smoothing  needs  to  be  developed 

IV  Coni  t  UNION 

The  frequency  domain  smoothing  approach  suggested  in  this 
correspondence  can  outperform  the  spatial  smoothing  approach. 


Angle  (degree) 

Fig  2  Beam  patterns  of  the  frequency  domain  smoothing  method  and  spa¬ 
tial  smoothing  method 

especially  for  small  sensor  arrays.  The  features  of  this  new  ap¬ 
proach  include  the  following: 

If  it  can  be  applied  to  arrays  with  arbitrary  geometry  and  sensor 
patterns ; 

2)  it  does  not  sacrifice  the  angular  resolution  for  smoothing;  and 

3)  it  docs  not  need  a  large  number  of  sensors  to  effect  the 
smoothing  operation. 

We  note  that  12|,  |7],  and  (8|  contain  some  interesting  results 
on  the  construction  of  the  approximations  to  the  linear  transform 
of  (3). 
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ABSTRACT 

This  paper  studies  the  effects  of  sensor  pat¬ 
tern  and  sensor  position  perturbations  on  the 
angle  estimation  performance  bound  for  multiple 
narrow-band  sources.  The  Cramer-Rao  Lower  Bound 
Is  used  with  a  probabilistic  modeling  of  the 
perturbations.  The  CRLBs  of  a  linear  uniform  ar¬ 
ray  under  sensor  position  and  pattern  perturba¬ 
tions  are  evaluated  In  detail  for  the  case  of  two 
narrow-band  sources. 

I.  INTRODUCTION 

In  many  applications  where  a  sensor  array  is 
used  for  angle  estimation  of  multiple  sources/ 
targets,  one  has  to  consider  the  effects  of  the 
sensor  pattern  and  sensor  position  perturbations. 
Though  some  efforts  have  been  made  to  develop 
robust  hlgh-resolution  estimates  [1|-(3|,  there  is 
a  need  to  understand  the  behavior  of  the  estima¬ 
tion  performance  bound  under  the  perturbations,  so 
that  one  can  see  how  much  room  left  for  developing 
more  robust  estimate. 

In  this  paper,  the  behavior  of  estimation 
performance  bound  under  sensor  position  and  pat¬ 
tern  perturbations  is  investigated.  The  Cramer- 
Rao  Lower  Bound  (CRLB)  is  chosen  since  it  is 
simple  to  evaluate,  and  since  all  available  high- 
resolutlon  direction  finding  methods  require  ei¬ 
ther  large  number  of  snapshots  or  high  SNR,  for 
which  tl«e  CRLB  has  been  seen  to  provide  a  rea¬ 
sonably  tight  bound  on  the  mean-square-error  (MSE) 
of  the  angle  estimate.  The  perturbations  on 
sensor  positions  and  patterns  are  modeled  as 
Gaussian  random  variables.  Ve  note  that  previous 
work  on  evaluation  of  the  CRLB  under  perturba¬ 
tions,  such  as  (4|-(6|,  involves  the  cases  of 
single  source  and  multiple  sources  disjoint  in  the 
frequency  domain. 

II.  MODEL  FORMULATION 

Consider  an  array  with  arbitrary  but  known 
nominal  geometry  as  shown  In  Fig.l.  Let  the  num¬ 
ber  of  sensors  be  M,  each  at  position  y1), 

y,) . ?„■(*„  y„),  with  respect  to 

the  chosen  reference  point.  Denote  the  sensor 
pattern  of  the  m-th  sensor  as  g„(9).  Let 

Ibis  work  was  supported  hy  RAOC  under  rtfr  Force 
Contract  F30602 -R0 -D-002 7 ( A-8- 1 1 24 ; 


? . be  the  wave-number  vectors  of  d 

narrow-band  plane  waves  of  the  same  central  fre¬ 
quency  f0.  The  array  output  complex  vector  X, 
Mx  1 ,  can  be  expressed  as 

X  .  AS  ♦  N  (  1  ) 

where  A  .  |a,  a,  . . . a,, I  Is  the  direction  matrix 

with  the  direction  vector 

at  -  (g1(6l)exp(J2n^l(-^1)) , 
gJ(6l)exp(J2rtf,(-'?1)) , 

gaOVexpUZniU-?,))!*,  1*1.2 . d,  (2) 

S  •  |S,  Sj  ...  S„|T  the  arriving  signal  vector 
modeled  as  zero-mean  Gaussian  random  vector  with  a 
diagonal  covariance  matrix,  and  N  -  |N,  N,  ...  N„|T 
the  receiver  noise  modeled  as  zero-mean  Gaussian 
random  vector  with  its  covariance  matrix  equal  to 
o*I.  S  and  N  are  assumed  to  be  independent. 

From  Fig.l,  we  have 

?„(-?,)  *  (x.cosD,  *  y»^ln9, ) / X  , 

m.1,2 . M,  1  =  1,2 . d.  (3) 

In  the  presence  of  sensor  position  per  turbo 1 1  on , 
the  m-th  sensor  position  vector  becomes 
jtm- ( x<*4x.  y„.6y„) ,  where  &xm  and  flym  are  the 
sensor  position  perturbations  In  x  direction  and  y 
direction  respectively.  The  sensor  position 
perturbations  6xm,  6y„,  m*l,2,...,M  ace  modeled  as 

1.1. d.  Gaussian  random  variables  of  zero  mean.  Ue 
assume  that  the  standard  deviation  of  the  position 
perturbations  is  much  smaller  than  the  spaces 
among  the  sensors. 

The  perturbed  complex  gains,  i.e.,  the 
perturbed  sensor  patterns,  are  modeled  as 
g.( 9) * ag.( 9)  where  g„(9),  m=l,2,...M  are  the 
nominal  patterns  and  flg„(9),  m»l,2,...H,  are 

1.1. d.  zero  mean  complex  Gaussian  random  perturba¬ 
tions.  Again,  the  standard  deviation  is  assumed 
t=  he  much  smaller  than  t lie  nominal  complex  gain 

R„  >■ 

In  order  to  see  the  amplitude  and  phase 
perturbation  effects  separately,  a  model  for  the 
case  of  pattern  amplitude  perturbation  onlv  1= 
also  set  up,  In  which  the  nominal  pattern  g,(  0)  , 
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a. 1,2 . H  are  assumed  to  be  positive  and  the 

perturbations  dig  „  (  0 ) ,  m-l,2....,M  to  be  1.1. d  real 
zero  mean  Caussian  random  variables. 

When  multiple  snapshots  are  used  for  estima¬ 
tion,  the  Independence  of  snapshots  is  assumed  and 
the  perturbations  are  assumed  to  remain  unchanged 
among  all  snapshots. 


where  90  is  the  parameter  set  at  the  nominal 
value.  For  Caussian  signals  and  noises  Ei| .  (10) 
becomes  |4] 


,  3R(0)  ,  3R(0) 

IJ.O)),  ,  =  trlR-‘(0)^-R-‘(0)-^r- 


e-90  (in 


where  R  Is  the  covariance  matrix  of  the  data  vec¬ 
tor  X. 


III.  CRAMER-RAO  LOWER  BOUND  UNDER  PERTURBATIONS 

If  we  only  consider  position  perturbations, 
l.e.,  assume  no  pattern  perturbations,  then  the 
unknown  real  parameter  set  9  contains 


9  -  (9,  e:  ...  9d  o‘, 
*x,  dy,  Ax,  Ay, 


9<,  »ii  °l, 


**„  *y„> 


l.e.,  a  total  of  2d»l»2M  real  parameters.  For  the 
case  of  pattern  perturbations  only,  the  unknown 
real  parameter  set  9  Is 

9  -  (0,  0;  ...  9,J  *>',  ... 

Re(  dgl(9l ) )  Imfdg/Q,))  ...  Re(Ag,(0d))  Im(Ag,(9d)) 
ReC  ag2ce, ) )  Im(  dg }( 9, ) )  ...  Re(Ag,(9d))  Im(Ag,(0d)) 

Re(*g„l9,))  I"'(  d«M(0, ) )  ...  KefAg^OJ)  I«.(dg„(Od))) 

(5) 

a  total  of  2d*l*2Md  real  parameters.  For  the  case 
of  pattern  amplitude  perturbations  only,  dgll(0), 
in « 1 , 2  ,  .  .  .  ,  H  are  set  up  as  1.1. d.  real  zero  mean 
Caussian  random  variables  In  section  II.  There¬ 
fore  the  unknown  real  parameter  set  In  this  case 
is 

0  .  (0,  0,  ...  ed  oj,  ...  o*d  o* 

*8,(9.)  *8,(9,)  *8,(9,) 

*8,(9,)  *8,(9,)  ...  *g,(9d) 


Because  the  first  2d*l  parameters  are  non- 
random,  J,(0)  has  it’s  first  2d*l  rows  and  columns 
equal  to  zero.  Let  Ap(JM<JI1|  be  the  covariance 
matrix  of  the  position  perturbations, 
the  covariance  matrix  of  the  pattern  perturbations 
(real  and  imaginary  part),  and  A.(nd,ddl  the 
covariance  matrix  of  the  pattern  amplitude  pertur¬ 
bation.  For  sensor  position  perturbation  only,  we 
have 


and  for  pattern  perturbation  only,  we  have 


J,(0)  . 


id.l  >  2d* l 


*9  I ,Md«  2nd  I 


and  for  pattern  amplitude  perturbation  only, 


■'2d.l  «  2d»l 


*g„(9,)  *«„(  9,) 


*8r,(9a>> 


*«IMd»ndl 


In  the  case  of  multiple  snapshots,  the  prob¬ 
ability  density  function  Pr(X|9)  In  Eq.(B)  would 
be  joint  density  function  Pr(  X, ,  X, ,  .  .  .  X„  |9),  where  K 
is  the  numbec  of  snapshots.  Under  the  assumption 
that  all  snapshots  are  Independently  identically 
distributed,  Pr ( X , , X , , . . . X„ |0)  .  [rr(X,|0)|K  and 
Eq.(8)  becomes 


a  total  of  2d»l*Hd  real  parameters. 

With  the  presence  of  random  parameters,  the 
Fisher  Information  matrix  Is  given  by  (7) 

J(9)  .  J , ( 9)  ♦  J,(9)  (7) 

where  the  elements  of  Jj  and  J2  are 

3Hog(Pr(X|0)) 

- >  <■> 


t  J2<s>  I  r.j  -  'V 


3Jlog(Pr(0) ) 


Under  the  assumption  that  the  position  and  pattern 
perturbations  are  small,  J,(9)  can  be  approximated 
by  the  corresponding  Fisher  information  matrix  for 
non-random  0  at  the  nominal  values  (4|,  l.e., 


IJ.tO)!,.,  =  -E,< 


3,log( Pr(X [9) ) 


|J.(0)I, .,  -  -K-E,  9( 


a’logtPrtX, | 9) ) 


The  J,(9)  will  remain  the  same  as  In  Eq.(9)  be¬ 
cause  of  the  assumption  that  the  perturbations 
remain  unchanged  among  all  snapshots. 

In  the  following  section,  we  will  numerically 
evaluate  the  CRLB  for  a  linear  uniform  array  to 
see  how  the  sensor  position  and  pattern  perturba¬ 
tions  affect  the  estimation  bound. 

IV.  CASE  STUDY 

Ve  consider  an  array  of  8  sensors  with  half 
wavelength  spacing,  which  presents  a  Rayleigh 
angular  resolution  of  16.4*  (degree).  For 
simplicity,  the  nominal  sensor  patterns  are  as¬ 
sumed  to  be  omnidirectional,  l.e.,  g(9).l  for  all 
sensors.  Let  9,  and  0,  be  the  angles  of  arrival 
of  two  narrow  band  sources  of  the  some  central 
frequency.  Again,  for  simplicity,  we  fix  9,  at 
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0*,  and  then  60  -  0,-0,  -  9,-  In  the  following 
figures,  the  CRL8  for  0,  estimation  under  pattern 
amplitude,  pattern  amplitude  and  phase,  and  posi¬ 
tion  perturbations,  denoted  as  CRLBa,  CRLBg  and 
CRLBp  respectively,  are  compared  vl th  the  CRL8  for 
0,  estimation  without  perturbations  (CRL8„).  For 
pattern  amplitude  perturbations  the  standard 
deviation  is  denoted  as  For  pattern  amplitude 
and  phase  perturbations,  which  are  modeled  as 
i.i.d.  zero  mean  complex  Gaussian  random  vari¬ 
ables,  we  denote  the  standard  deviation  of  its 
real  part  or  imaginary  part  as  o?,i.e.,the  stan¬ 
dard  deviation  of  the  complex  pattern  perturbation 
is  For  sensor  position  perturbations,  we 
let  the  standard  deviation  of  the  perturbations  in 
x  direction  and  y  direction  be  the  same  and 
denoted  as  o  . 

P 

Figure  2  shows  the  effect  of  the  pattern 
amplitude  perturbations  on  the  bound  (CRLBa)  as  a 
function  of  00,  the  angular  separation  between  the 
two  narrow-band  sources.  For  00  larger  than  the 
resolution  cell  (=16.4'),  the  CRLBa  is  almost  the 
same  as  CRLB0,  i.e.,  the  amplitude  perturbation  of 
the  sensor  gains  has  almost  no  effect  on  the 
estimation  performance  potential.  In  the  high- 
resolution  region  (09<16.4*),  however,  the  larger 
the  amplitude  perturbation,  the  poorer  the  estima¬ 
tion  performance  potential,  which  implies  some  ex¬ 
tra  difficulty  to  overcome  in  order  to  achieve 
high  resolution. 

To  overcome  such  an  extra  difficulty  for  high 
resolution,  Fig.  3  shows  that  one  possible  way  is 
to  increase  the  signal-to-noise  ratio  (SNR).  In 
this  sense  we  may  loosely  consider  that  the  pat¬ 
tern  amplitude  perturbation  has  a  similar  effect 
on  the  performance  potential  as  the  receiver 
noise. 

Figure  4,  Figure  5  and  Figure  6  shov  the  CRLB 
under  both  the  amplitude  and  phase  perturbation  on 
the  sensor  pattern  (CRLBg).  When  perturbations  on 
the  phases  of  the  complex  sensor  gain  are  added, 
the  CRLBg  behaves  quite  differently  from  that 
without  phase  perturbation.  Figure  4  shows  the 
bound  as  a  function  of  the  angle  separation  00 
with  different  e  ,  while  Figure  5  as  of  SNR  with 
different  00  and  Figure  6  as  of  SNR  with  different 
<j?.  Both  Figure  5  and  Figure  6  indicate  that  with 
the  phase  perturbations  on  the  sensor  gain,  the 
CRLBg  will  level  off  as  SNR  increases.  That  is, 
the  CRLBg  can  not  be  reduced  to  arbitrarily  small 
by  increasing  SNR.  Under  sensor  pattern  perturba¬ 
tion,  therefore,  it  is  the  phase  perturbation  that 
significantly  affects  the  estimation  performance. 

Figures  2-9  show  the  CRLB  under  sensor  posi¬ 
tion  perturbations  (CRLBp).  We  can  see  that  the 
effects  of  position  perturbation  are  almost  the 
same  as  that  of  pattern  perturbation.  Again,  when 
noise  is  no  longer  a  dominant  factor,  it  is  the 
position  perturbation  that  will  limit  the  estima¬ 
tion  performance. 


Fig.  7  plots  CRLBp  as  a  function  of  the 
angular  separation  09.  In  the  1  ow- r eso i u t  i  on 
region  (68>16.4')t  the  degradation  from  the  CRLB, 
is  almost  independent  of  00  and  determined  only  by 
the  position  perturbation  standard  deviation  op . 
In  the  high-resolution  region  (60C16.4'),  however, 
the  degradation  from  CRLB,  is  00-dependent  and  be¬ 
comes  smaller  when  00  is  smaller.  Noting  the 
sharp  Increase  of  CRLB,  in  the  high-resolution 
region,  we  should  realize  that  the  smaller 
degradation  with  smaller  00  merely  means  t lie  mote 
dominant  factor  which  the  receiver  noise  shows 
with  smaller  00. 

Fig.  10  shows  the  effect  of  the  number  of 
snapshots  on  the  CRLB  under  position  perturbations 
(CRLBp),  from  which  the  Increase  of  the  number  of 
snapshots  is  seen  to  be  similar  to  the  increase  of 
the  signal  to  noise  ratio.  That  is,  <ns  the  number 
of  snapshots  increases  the  CRLBp  levels  off. 

V.  CONCLUSIONS 

Both  sensor  pattern  and  position  perturbation 
can  seriously  degrade  the  potential  estimation 
performance.  With  proper  joint  estimation  scheme, 
the  effect  of  the  pattern  amplitude  perturbation 
might  be  compensated  by  Increasing  the  signal  to 
noise  ratio.  When  the  signal  to  noise  ratio  is 
high  the  position  perturbation  or  pattern  phase 
perturbation  becomes  the  dominant  factor  and  the 
CRLBp  or  CRLBg  will  level  off.  In  such  situa¬ 
tions,  one  might  not  be  able  to  reduce  the  CRLng 
and  CRLBp  to  an  arbitrarily  small  number  by  fur¬ 
ther  increase  of  the  signal  to  noise  ratio,  even 
with  some  joint  estimation  scheme.  Under  the  pat¬ 
tern  or  position  perturbation,  the  increase  of  the 
number  of  snapshots  is  expected  to  offer  only 
limited  help  as  the  SNR. 
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